Vibrio-Sequins - dPCR-traceable DNA standards for quantitative genomics of Vibrio spp.

BACKGROUND: Vibrio spp. are a diverse group of ecologically important marine bacteria responsible for several foodborne outbreaks of gastroenteritis around the world. Their detection and characterization are moving away from conventional culture-based methods towards next generation sequencing (NGS)-based approaches. However, genomic methods are relative in nature and suffer from technical biases arising from library preparation and sequencing. Here, we introduce a quantitative NGS-based method that enables the quantitation of Vibrio spp. at the limit of quantification (LOQ) through artificial DNA standards and their absolute quantification via digital PCR (dPCR). RESULTS: We developed six DNA standards, called Vibrio-Sequins, together with optimized TaqMan assays for their quantification in individually sequenced DNA libraries via dPCR. To enable Vibrio-Sequin quantification, we validated three duplex dPCR methods to quantify the six targets. LOQs were ranging from 20 to 120 cp/µl for the six standards, whereas the limit of detection (LOD) was ~ 10 cp/µl for all six assays. Subsequently, a quantitative genomics approach was applied to quantify Vibrio-DNA in a pooled DNA mixture derived from several Vibrio species in a proof-of-concept study, demonstrating the increased power of our quantitative genomic pipeline through the coupling of NGS and dPCR. CONCLUSIONS: We significantly advance existing quantitative (meta)genomic methods by ensuring metrological traceability of NGS-based DNA quantification. Our method represents a useful tool for future metagenomic studies aiming at quantifying microbial DNA in an absolute manner. The inclusion of dPCR into sequencing-based methods supports the development of statistical approaches for the estimation of measurement uncertainties (MU) for NGS, which is still in its infancy.

Adaptive evolution and segregating load contribute to the genomic landscape of divergence in two tree species connected by episodic gene flow.

Speciation often involves repeated episodes of genetic contact between divergent populations before reproductive isolation (RI) is complete. Whole-genome sequencing (WGS) holds great promise for unravelling the genomic bases of speciation. We have studied two ecologically divergent, hybridizing species of the 'model tree' genus Populus (poplars, aspens, cottonwoods), Populus alba and P. tremula, using >8.6 million single nucleotide polymorphisms (SNPs) from WGS of population pools. We used the genomic data to (i) scan these species' genomes for regions of elevated and reduced divergence, (ii) assess key aspects of their joint demographic history based on genomewide site frequency spectra (SFS) and (iii) infer the potential roles of adaptive and deleterious coding mutations in shaping the genomic landscape of divergence. We identified numerous small, unevenly distributed genome regions without fixed polymorphisms despite high overall genomic differentiation. The joint SFS was best explained by ancient and repeated gene flow and allowed pinpointing candidate interspecific migrant tracts. The direction of selection (DoS) differed between genes in putative migrant tracts and the remainder of the genome, thus indicating the potential roles of adaptive divergence and segregating deleterious mutations on the evolution and breakdown of RI. Genes affected by positive selection during divergence were enriched for several functionally interesting groups, including well-known candidate 'speciation genes' involved in plant innate immunity. Our results suggest that adaptive divergence affects RI in these hybridizing species mainly through intrinsic and demographic processes. Integrating genomic with molecular data holds great promise for revealing the effects of particular genetic pathways on speciation.

Selection against recombinant hybrids maintains reproductive isolation in hybridizing Populus species despite F1 fertility and recurrent gene flow.

Natural hybrid zones have proven to be precious tools for understanding the origin and maintenance of reproductive isolation (RI) and therefore species. Most available genomic studies of hybrid zones using whole- or partial-genome resequencing approaches have focused on comparisons of the parental source populations involved in genome admixture, rather than exploring fine-scale patterns of chromosomal ancestry across the full admixture gradient present between hybridizing species. We have studied three well-known European 'replicate' hybrid zones of Populus alba and P. tremula, two widespread, ecologically divergent forest trees, using up to 432 505 single-nucleotide polymorphisms (SNPs) from restriction site-associated DNA (RAD) sequencing. Estimates of fine-scale chromosomal ancestry, genomic divergence and differentiation across all 19 poplar chromosomes revealed strikingly contrasting results, including an unexpected preponderance of F1 hybrids in the centre of genomic clines on the one hand, and genomically localized, spatially variable shared variants consistent with ancient introgression between the parental species on the other. Genetic ancestry had a significant effect on survivorship of hybrid seedlings in a common garden trial, pointing to selection against early-generation recombinants. Our results indicate a role for selection against recombinant genotypes in maintaining RI in the face of apparent F1 fertility, consistent with the intragenomic 'coadaptation' model of barriers to introgression upon secondary contact. Whole-genome resequencing of hybridizing populations will clarify the roles of specific genetic pathways in RI between these model forest trees and may reveal which loci are affected most strongly by its cyclic breakdown.

Genome-wide patterns of differentiation and spatially varying selection between postglacial recolonization lineages of Populus alba (Salicaceae), a widespread forest tree.

Studying the divergence continuum in plants is relevant to fundamental and applied biology because of the potential to reveal functionally important genetic variation. In this context, whole-genome sequencing (WGS) provides the necessary rigour for uncovering footprints of selection. We resequenced populations of two divergent phylogeographic lineages of Populus alba (n = 48), thoroughly characterized by microsatellites (n = 317), and scanned their genomes for regions of unusually high allelic differentiation and reduced diversity using > 1.7 million single nucleotide polymorphisms (SNPs) from WGS. Results were confirmed by Sanger sequencing. On average, 9134 high-differentiation (≥ 4 standard deviations) outlier SNPs were uncovered between populations, 848 of which were shared by ≥ three replicate comparisons. Annotation revealed that 545 of these were located in 437 predicted genes. Twelve percent of differentiation outlier genome regions exhibited significantly reduced genetic diversity. Gene ontology (GO) searches were successful for 327 high-differentiation genes, and these were enriched for 63 GO terms. Our results provide a snapshot of the roles of 'hard selective sweeps' vs divergent selection of standing genetic variation in distinct postglacial recolonization lineages of P. alba. Thus, this study adds to our understanding of the mechanisms responsible for the origin of functionally relevant variation in temperate trees.

Generation of SNP datasets for orangutan population genomics using improved reduced-representation sequencing and direct comparisons of SNP calling algorithms.

BACKGROUND: High-throughput sequencing has opened up exciting possibilities in population and conservation genetics by enabling the assessment of genetic variation at genome-wide scales. One approach to reduce genome complexity, i.e. investigating only parts of the genome, is reduced-representation library (RRL) sequencing. Like similar approaches, RRL sequencing reduces ascertainment bias due to simultaneous discovery and genotyping of single-nucleotide polymorphisms (SNPs) and does not require reference genomes. Yet, generating such datasets remains challenging due to laboratory and bioinformatical issues. In the laboratory, current protocols require improvements with regards to sequencing homologous fragments to reduce the number of missing genotypes. From the bioinformatical perspective, the reliance of most studies on a single SNP caller disregards the possibility that different algorithms may produce disparate SNP datasets. RESULTS: We present an improved RRL (iRRL) protocol that maximizes the generation of homologous DNA sequences, thus achieving improved genotyping-by-sequencing efficiency. Our modifications facilitate generation of single-sample libraries, enabling individual genotype assignments instead of pooled-sample analysis. We sequenced ~1% of the orangutan genome with 41-fold median coverage in 31 wild-born individuals from two populations. SNPs and genotypes were called using three different algorithms. We obtained substantially different SNP datasets depending on the SNP caller. Genotype validations revealed that the Unified Genotyper of the Genome Analysis Toolkit and SAMtools performed significantly better than a caller from CLC Genomics Workbench (CLC). Of all conflicting genotype calls, CLC was only correct in 17% of the cases. Furthermore, conflicting genotypes between two algorithms showed a systematic bias in that one caller almost exclusively assigned heterozygotes, while the other one almost exclusively assigned homozygotes. CONCLUSIONS: Our enhanced iRRL approach greatly facilitates genotyping-by-sequencing and thus direct estimates of allele frequencies. Our direct comparison of three commonly used SNP callers emphasizes the need to question the accuracy of SNP and genotype calling, as we obtained considerably different SNP datasets depending on caller algorithms, sequencing depths and filtering criteria. These differences affected scans for signatures of natural selection, but will also exert undue influences on demographic inferences. This study presents the first effort to generate a population genomic dataset for wild-born orangutans with known population provenance.

Genomic scan for single nucleotide polymorphisms reveals patterns of divergence and gene flow between ecologically divergent species.

Recent advances in population genomics have triggered great interest in the genomic landscape of divergence in taxa with 'porous' species boundaries. One important obstable of previous studies of this topic was the low genomic coverage achieved. This issue can now be overcome by the use of 'next generation' or short-read DNA-sequencing approaches capable of assaying many thousands of single nucleotide polymorphisms (SNPs) in divergent species. We have scanned the 'porous' genomes of Populus alba and Populus tremula, two ecologically divergent hybridizing forest trees, using >38,000 SNPs assayed by restriction site associated DNA (RAD) sequencing. Windowed analyses indicate great variation in genetic divergence (e.g. the proportion of fixed SNPs) between species, and these results are unlikely to be strongly biased by genomic features of the Populus trichocarpa reference genome used for SNP calling. Divergence estimates were significantly autocorrelated (P < 0.01; Moran's I up to 0.6) along 11 of 19 chromosomes. Many of these autocorrelations involved low divergence blocks, thus suggesting that allele sharing was caused by recurrent gene flow rather than shared ancestral polymorphism. A conspicuous low divergence block of three megabases was detected on chromosome XIX, recently put forward as an incipient sex chromosome in Populus, and was largely congruent with introgression of mapped microsatellites in two natural hybrid zones (N > 400). Our results help explain the origin of the 'genomic mosaic' seen in these taxa with 'porous' genomes and suggest rampant introgression or extensive among-species conservation of an incipient plant sex chromosome. RAD sequencing holds great promise for detecting patterns of divergence and gene flow in highly divergent hybridizing species.

Cost-effective fluorescent amplified fragment length polymorphism (AFLP) analyses using a three primer system.

The amplified fragment length polymorphism (AFLP) technique is a widely used multi-purpose DNA fingerprinting tool. The ability to size-separate fluorescently labelled AFLP fragments on a capillary electrophoresis instrument has provided a means for high-throughput genome screening, an approach particularly useful in studying the molecular ecology of nonmodel organisms. While the 'per-marker-generated' costs for AFLP are low, fluorescently labelled oligonucleotides remain costly. We present a cost-effective method for fluorescently end-labelling AFLPs that should make this tool more readily accessible for laboratories with limited budgets. Both standard fluorescent AFLPs and the end-labelled alternatives presented here are repeatable and produce similar numbers of fragments when scored using both manual and automated scoring methods. While it is not recommended to combine data using the two approaches, the results of the methods are qualitatively comparable, indicating that AFLP end-labelling is a robust alternative to standard methods of AFLP genotyping. For researchers commencing a new AFLP project, the AFLP end-labelling method outlined here is easily implemented, as it does not require major changes to PCR protocols and can significantly reduce the costs of AFLP studies.

Eukaryotic transcriptomics in silico: optimizing cDNA-AFLP efficiency.

BACKGROUND: Complementary-DNA based amplified fragment length polymorphism (cDNA-AFLP) is a commonly used tool for assessing the genetic regulation of traits through the correlation of trait expression with cDNA expression profiles. In spite of the frequent application of this method, studies on the optimization of the cDNA-AFLP assay design are rare and have typically been taxonomically restricted. Here, we model cDNA-AFLPs on all 92 eukaryotic species for which cDNA pools are currently available, using all combinations of eight restriction enzymes standard in cDNA-AFLP screens. RESULTS: In silco simulations reveal that cDNA pool coverage is largely determined by the choice of individual restriction enzymes and that, through the choice of optimal enzyme combinations, coverage can be increased from <40% to 75% without changing the underlying experimental design. We find evidence of phylogenetic signal in the coverage data, which is largely mediated by organismal GC content. There is nonetheless a high degree of consistency in cDNA pool coverage for particular enzyme combinations, indicating that our recommendations should be applicable to most eukaryotic systems. We also explore the relationship between the average observed fragment number per selective AFLP-PCR reaction and the size of the underlying cDNA pool, and show how AFLP experiments can be used to estimate the number of genes expressed in a target tissue. CONCLUSION: The insights gained from in silico screening of cDNA-AFLPs from a broad sampling of eukaryotes provide a set of guidelines that should help to substantially increase the efficiency of future cDNA-AFLP experiments in eukaryotes. In silico simulations also suggest a novel use of cDNA-AFLP screens to determine the number of transcripts expressed in a target tissue, an application that should be invaluable as next-generation sequencing technologies are adapted for differential display.

Distinct migratory and non-migratory ecotypes of an endemic New Zealand eleotrid (Gobiomorphus cotidianus) - implications for incipient speciation in island freshwater fish species.

BACKGROUND: Many postglacial lakes contain fish species with distinct ecomorphs. Similar evolutionary scenarios might be acting on evolutionarily young fish communities in lakes of remote islands. One process that drives diversification in island freshwater fish species is the colonization of depauperate freshwater environments by diadromous (migratory) taxa, which secondarily lose their migratory behaviour. The loss of migration limits dispersal and gene flow between distant populations, and, therefore, is expected to facilitate local morphological and genetic differentiation. To date, most studies have focused on interspecific relationships among migratory species and their non-migratory sister taxa. We hypothesize that the loss of migration facilitates intraspecific morphological, behavioural, and genetic differentiation between migratory and non-migratory populations of facultatively diadromous taxa, and, hence, incipient speciation of island freshwater fish species. RESULTS: Microchemical analyses of otolith isotopes (88Sr, 137Ba and 43Ca) differentiated migratory and non-migratory stocks of the New Zealand endemic Gobiomorphus cotidianus McDowall (Eleotridae). Samples were taken from two rivers, one lake and two geographically-separated outgroup locations. Meristic analyses of oculoscapular lateral line canals documented a gradual reduction of these structures in the non-migratory populations. Amplified fragment length polymorphism (AFLP) fingerprints revealed considerable genetic isolation between migratory and non-migratory populations. Temporal differences in reproductive timing (migratory = winter spawners, non-migratory = summer spawners; as inferred from gonadosomatic indices) provide a prezygotic reproductive isolation mechanism between the two ecotypes. CONCLUSION: This study provides a holistic look at the role of diadromy in incipient speciation of island freshwater fish species. All four analytical approaches (otolith microchemistry, morphology, spawning timing, population genetics) yield congruent results, and provide clear and independent evidence for the existence of distinct migratory and non-migratory ecotypes within a river in a geographically confined range. The morphological changes within the non-migratory populations parallel interspecific patterns observed in all non-migratory New Zealand endemic Gobiomorphus species and other derived gobiid taxa, a pattern suggesting parallel evolution. This study indicates, for the first time, that distinct ecotypes of island freshwater fish species may be formed as a consequence of loss of migration and subsequent diversification. Therefore, if reproductive isolation persists, these processes may provide a mechanism to facilitate speciation.

Mate choice for genetic quality when environments vary: suggestions for empirical progress.

Mate choice for good-genes remains one of the most controversial evolutionary processes ever proposed. This is partly because strong directional choice should theoretically deplete the genetic variation that explains the evolution of this type of female mating preference (the so-called lek paradox). Moreover, good-genes benefits are generally assumed to be too small to outweigh opposing direct selection on females. Here, we review recent progress in the study of mate choice for genetic quality, focussing particularly on the potential for genotype by environment interactions (GEIs) to rescue additive genetic variation for quality, and thereby resolve the lek paradox. We raise five questions that we think will stimulate empirical progress in this field, and suggest directions for research in each area: (1) How is condition-dependence affected by environmental variation? (2) How important are GEIs for maintaining additive genetic variance in condition? (3) How much do GEIs reduce the signalling value of male condition? (4) How does GEI affect the multivariate version of the lek paradox? (5) Have mating biases for high-condition males evolved because of indirect benefits?

Male pregnancy in seahorses and pipefish: beyond the mammalian model.

Pregnancy has been traditionally defined as the period during which developing embryos are incubated in the body after egg-sperm union. Despite strong similarities between viviparity in mammals and other vertebrate groups, researchers have historically been reluctant to use the term pregnancy for non-mammals in recognition of the highly developed form of viviparity in eutherians. Syngnathid fishes (seahorses and pipefishes) have a unique reproductive system, where the male incubates developing embryos in a specialized brooding structure in which they are aerated, osmoregulated, protected and likely provisioned during their development. Recent insights into physiological, morphological and genetic changes associated with syngnathid reproduction provide compelling evidence that male incubation in these species is a highly specialized form of reproduction akin to other forms of viviparity. Here, we review these recent advances, highlighting similarities and differences between seahorse and mammalian pregnancy. Understanding the changes associated with the parallel evolution of male pregnancy in the two major syngnathid lineages will help to identify key innovations that facilitated the development of this unique form of reproduction and, through comparison with other forms of live bearing, may allow the identification of a common set of characteristics shared by all viviparous organisms.

Monitoring the effects of pulp and paper effluent is restricted in genetically distinct populations of common bully (Gobiomorphus cotidianus).

The common bully (Gobiomorphus cotidianus), a small-bodied New Zealand native fish species, was used to monitor population impacts of multiple effluents in the Tarawera River, New Zealand. In an initial survey, the absence of reproductive development at the expected spawning time for common bullywas observed in a population downstream of effluent discharges. Subsequently, we examined the hypotheses that the observed changes were due to effluent exposure, migratory patterns, or genetic differences between populations. Liver detoxification enzyme activity and stable isotopes provided evidence against upstream migration of sexually mature bully. The observed presence of developed gonads in the downstream population during winter season resulted in the rejection of the hypothesis that reproductive failure was due to effluent exposure, and itwas concluded that there were substantial differences in reproductive timing. Genetic analyses of two upstream, one downstream, and one population from a nearby coastal river indicated the upstream (reference) and downstream (effluent exposed) bully in the river formed genetically distinct populations. The identification of a nearby river population with similar reproductive timing and high genetic similarity to the effluent-exposed population suggests that the observed differences in the genetics of the downstream population were not caused by effluent exposure. The genetic analysis did highlight the lack of downstream dispersion and gene flow in the river which could possibly be related to anthropogenic stress.

A BAC library for the goldfish Carassius auratus auratus (Cyprinidae, Cypriniformes).

A goldfish (Carassius auratus auratus) bacterial artificial chromosome genomic library (BAC library) was constructed from one aquarium-bred male specimen (tetraploid, 4n=100, genome size=3.52 pg/cell). The library consists of 128,352 positive clones with an average insert size of 150.4 kb, covering the genome 11-fold. All clones were spotted onto nylon filters and thus are available for screening of genomic regions of interest, such as candidate genes, gene families, or large-sized syntenic DNA regions of cyprinid species. Preliminary screens with two genes were conducted with hybridizing probes to the genes RAG1 and lgi1. RAG1 is a single-copy gene in zebrafish and is duplicated in C. a. auratus. We found a very close correlation between the number of positive BAC clones and the expected library coverage. Two copies of lgi1 were found in zebrafish. We have detected four different copies in C. a. auratus, not in the expected abundance, which indicates some variation in the coverage of the BAC library. The preliminary screens indicate that many duplicated genes that resulted from the ancient fish-specific genome duplication persist in the tetraploid goldfish genome. Hence, the BAC library will provide a useful resource for the future work on comparative genomics, polyploidy, diploidization, and evolutionary genomics in fishes.

Sympatric speciation in Nicaraguan crater lake cichlid fish.

Sympatric speciation, the formation of species in the absence of geographical barriers, remains one of the most contentious concepts in evolutionary biology. Although speciation under sympatric conditions seems theoretically possible, empirical studies are scarce and only a few credible examples of sympatric speciation exist. Here we present a convincing case of sympatric speciation in the Midas cichlid species complex (Amphilophus sp.) in a young and small volcanic crater lake in Nicaragua. Our study includes phylogeographic, population-genetic (based on mitochondrial DNA, microsatellites and amplified fragment length polymorphisms), morphometric and ecological analyses. We find, first, that crater Lake Apoyo was seeded only once by the ancestral high-bodied benthic species Amphilophus citrinellus, the most common cichlid species in the area; second, that a new elongated limnetic species (Amphilophus zaliosus) evolved in Lake Apoyo from the ancestral species (A. citrinellus) within less than approximately 10,000 yr; third, that the two species in Lake Apoyo are reproductively isolated; and fourth, that the two species are eco-morphologically distinct.

Plate tectonics and biogeographical patterns of the Pseudophoxinus (Pisces: Cypriniformes) species complex of central Anatolia, Turkey.

We investigated the phylogenetic relationships of Pseudophoxinus (Cyprinidae: Leuciscinae) species from central Anatolia, Turkey to test the hypothesis of geographic speciation driven by early Pliocene orogenic events. We analyzed 1141 aligned base pairs of the complete cytochrome b mitochondrial gene. Phylogenetic relationships reconstructed by maximum likelihood, Bayesian likelihood, and maximum parsimony methods are identical, and generally well supported. Species and clades are restricted to geologically well-defined units, and are deeply divergent from each other. The basal diversification of central Anatolian Pseudophoxinus is estimated to have occurred approximately 15 million years ago. Our results are in agreement with a previous study of the Anatolian fish genus Aphanius that also shows a diversification pattern driven by the Pliocene orogenic events. The distribution of clades of Aphanius and Pseudophoxinus overlap, and areas of distribution comprise the same geological units. The geological history of Anatolia is likely to have had a major impact on the diversification history of many taxa occupying central Anatolia; many of these taxa are likely to be still unrecognized as distinct.

Molecular phylogeny and historical biogeography of the Aphanius (Pisces, Cyprinodontiformes) species complex of central Anatolia, Turkey.

Phylogenetic relationships of a subset of Aphanius fish comprising central Anatolia, Turkey, are investigated to test the hypothesis of geographic speciation driven by early Pliocene orogenic events in spite of morphological similarity. We use 3434 aligned base pairs of mitochondrial DNA from 42 samples representing 36 populations of three species and six outgroup species to test this hypothesis. Genes analyzed include those encoding the 12S and 16S ribosomal RNAs; transfer RNAs coding for valine, leucine, isoleucine, glutamine, methionine, tryptophan, alanine, asparagine, cysteine, and tyrosine; and complete NADH dehydrogenase subunits I and II. Distance based minimum evolution and maximum-likelihood analyses identify six well-supported clades consisting of Aphanius danfordii, Aphanius sp. aff danfordii, and four clades of Aphanius anatoliae. Parsimony analysis results in 462 equally parsimonious trees, all of which contain the six well supported clades identified in the other analyses. Our phylogenetic results are supported by hybridization studies (Villwock, 1964), and by the geological history of Anatolia. Phylogenetic relationships among the six clades are only weakly supported, however, and differ among analytical methods. We therefore test and subsequently reject the hypothesis of simultaneous diversification among the six central Anatolian clades. However, our analyses do not identify any internodes that are significantly better supported than expected by chance alone. Therefore, although bifurcating branching order is hypothesized to underlie this radiation, the exact branching order is difficult to estimate with confidence.